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Abstract. Effective electron mobilities are obtained by transport measurements 

on InAs nanowire field-effect transistors at temperatures ranging from 10 — 200 K. 
The mobility-versus-temperature is found to have a positive slope at temperatures 
below ~ 30 — 50 K, which turns over to a negative slope at higher temperatures. 
The magnitude and temperature dependence of the observed mobility are both 
explained by the identification of Coulomb scattering from ionized surface states as 
the dominant scattering mechanism. The negative slope above 50 K is primarily due 
to the thermally activated increase in the number of scatterers, although nanoscale 
confinement also plays a role as higher radial subbands are populated, leading to 
interband scattering and a shift of the carrier distribution closer to the surface. 
Scattering rate calculations using finite-element simulations of the nanowire transistor 
confirm that these mechanisms can fully explain the experimental data. 

PACS numbers: 73.63.-b 
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1. Introduction 

Semiconductor nanowires grown by the vapour-liquid-solid (VLS) method [H El El SI IS] 
are the subject of active study, with many potential applications ranging from nanoscale 
circuits [6J and gas sensors [7J to high-efficiency solar cells [HI El [IHl E] • In particular, 
InAs nanowires form Ohmic contacts easily [T2], and can be grown with low structural 
defect densities [3] , giving rise to high electron mobilities [13] , though still low compared 
to high-quality bulk InAs [H]. The quasi-one- dimensional nature of electron transport 
at low temperatures [15] together with a spin-orbit coupling ~ 40 times larger than 
GaAs makes InAs an attractive material for the development of spintronic devices such 
as electron spin qubits in gate-defined quantum dots [Ml |T71 118|. Although transport 
in InAs nanowires is well-studied [I^, the detailed role played by surface states and 
the surface potential [20l EH [221 EH] with regard to the electron mobility is not well 
understood. 

In this paper, we present temperature dependent electron mobility measurements 
on InAs nanowire field-effect transistors (FETs) that show a characteristic temperature 
dependence across a wide range of nanowire diameters. The mobility peaks in the range 
3, 000 — 20, 000 cm^V^^s"^ near 40 K, with a positive slope at lower temperatures and 
a negative slope at higher temperatures. These mobilities are too low to be limited by 
optical or acoustic phonon scattering, which yield estimated bulk mobilities 2 — 3 orders 
of magnitude larger than those observed here [2l]. We expect this to remain true even 
in quasi-one-dimensional systems, where phonon scattering is moderately enhanced due 
to a larger available phase space for scattering [25]. The nanowire geometry suggests 
that surface scattering of some form should be dominant. Surface states are known to 
be present at densities ~ 10^^ — 10^^ cm^^ eV~^ and to act as donors. We argue 
that a sufficient density of these positively charged surface states should be more 
effective at scattering electrons than surface roughness, and therefore limit the mobility. 
Consistent with this picture, chemical treatment of the nanowire surface is seen to 
have a strong effect on the temperature-dependent mobility. To further confirm this 
hypothesis, we calculate Coulomb scattering rates using a Poisson-Schrodinger solver to 
carry out numerical simulations of the FET devices including surface charges. We find 
that the decrease in mobility with temperature above ~ 50 K can only arise from an 
increase in the number of ionized surface states, presumably due to thermal activation. 
Although populating higher radial subbands with increasing temperature contributes to 
an increase in the scattering rate, it is not alone sufficient to cause a negative slope in 
the mobility versus temperature. Our results underscore the need for tailored surface 
passivation techniques [261 EZ] to reduce the density of surface scatterers and smooth the 
local electronic potential, leading to increased carrier mobility and more ideal devices 
for a wide range of quantum transport and optoelectronics applications. These results 
on pure InAs nanowires serve as a benchmark to compare with studies of passivated 
nanowires. 
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Figure 1. (a) Low and (b) high magnification bright-field TEM images of an InAs 
nanowire grown by GS-MBE at 0.5 fim/hr. Scale bars are 500 nm in (a) and 5 nm 
in (b). The inset in (b) shows selected area diffraction pattern along the [2110] zone 
axis indicating pure wurtzite crystal structure. A majority of wires grown under these 
conditions had low stacking fault densities < 1 /im~^. 



2. Nanowire growth by gas-source MBE 

InAs nanowires were grown in a gas source molecular beam epitaxy (GS-MBE) system 
using Au seed particles. A 1 nm Au film is heated to form nanoparticles on a GaAs 
(lll)B substrate. For nanowire growth, In atoms were supplied as monomers from an 
effusion cell, and As2 dimers were supplied from an AsHs gas cracker operating at 950°C. 
Nanowire growth proceeded at a substrate temperature of 420°C, an In impingement 
rate of 0.5 /im/hr, and a V/III flux ratio of 4. The nanowires grew in random orientations 
with respect to the GaAs (lll)B substrate, possibly due to the large lattice mismatch 
strain between InAs and GaAs. Transmission electron microscopy (TEM) analysis, 
shown in [l^, indicated a Au nanoparticle at the end of each nanowire (darker contrast 
at the left end), consistent with the VLS process. Most nanowires had a rod-shaped 
morphology with negligible tapering and a diameter (~ 20 — 80 nm) that was roughly 
equal to the Au nanoparticle diameter at the top of each nanowire, indicating minimal 
sidewall deposition. 

A common occurrence in III-V nanowires is the existence of stacking faults whereby 
the crystal structure alternates between zincblende and wurtzite, or exhibits twinning, 
along the nanowire length. Joyce et al. |1] and Dick et al. [5] have shown that growth 
parameters in metalorganic chemical vapour deposition (MOCVD) have profound effects 
on the InAs nanowire crystal phase. Zincblende, wurtzite, or mixed zincblende/wurtzite 
nanowires were formed by simply tuning the temperature and V/III ratio. We have 
found that for GS-MBE grown InAs nanowires, stacking faults can be nearly eliminated 
and pure wurtzite structures can be realized at sufficiently low growth rate ~ 0.5/im/hr. 
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Figure 2. (a-c) Conductance versus backgate voltage for devices 1 — 3 at selected 
temperatures. D is the nanowire diameter and L the FET channel length (device 1 is 
tapered with an average nanowire diameter (I?) = 71 nm). The tangent lines drawn on 
the T = 122 K and T = 60 K traces in (a) indicate the maximum slopes corresponding 
to peak field-effect mobility. The pinchofF threshold voltage is defined as the intercept 
between this tangent line and the G = axis, (d-f) The pinchofF threshold voltages 
versus temperature extracted from the conductance measurements. In (d), data are 
shown for device 1 before and after an ammonium sulfide treatment was applied to the 
FET channel (the data in (a) correspond to the untreated case). The empirical fits in 
(d-f) are of the form Vt = Vq + VxeT^^I^'^ ^ as described in the text. 



At higher growth rates, but otherwise identical growth conditions, the InAs nanowires 
exhibited a much larger fraction of stacking faults on average. For example, TEM 
analysis of InAs nanowires grown at a rate of 1 /im hr~^ exhibited an average linear 
density of stacking faults ~ 1 /im"^. Similar to GaAs nanowires pHl |2Hl S], the density 
of faults diminished dramatically when the growth rate was reduced. Selected area 
electron diffraction for a typical nanowire (inset of [Tjo) confirms the pure wurtzite crystal 
structure and the absence of stacking faults. 
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3. Mobility in field-effect transistors 

Field-effect transistors (FETs) were fabricated by meclianically depositing as-grown 
nanowires on a 175 nm tliick Si02 layer above a n^-Si substrate that functions as 
a backgate, and writing source/drain contacts for selected wires using electron-beam 
lithography. This was followed by an ammonium sulfide etching and chemical passiva- 
tion process to remove the native oxide and prevent regrowth [12] prior to evaporation 
of Ni/Au contacts. Channel lengths ranged from 0.7-3 fim. Transport measurements 
were carried out in He vapour in an Oxford continuous flow cryostat from 4 K to room 
temperature. Bias and gate voltages were applied using a high resolution home-built 
voltage source, and a DL Instruments current preamplifier was used to measure DC cur- 
rent at a noise floor ~ 0.5 pA/\/ Hz. All devices tested at room temperature displayed 
fully Ohmic I-V characteristics, with resistances typically in the range of 5 — 200 kfl. 



3.1. Results 

We investigated many such devices to varying levels of detail, and found qualitatively 
similar results. Here we will focus on three representative devices, denoted 1, 2 and 3 
with nanowire diameters D = 71, 50 and 35 nm, respectively. The nanowires in devices 

2 and 3 were untapered, whereas the nanowire in device 1 was tapered, with diameter 
linearly varying from 53 nm to 90 nm across the FET channel (average diameter {D) = 
71 nm). Subsequent to transport measurements, the channel of device 1 was subjected 
to an ammonium sulfide etching and passivation treatment similar to that carried out 
just before contacting. Its temperature dependent mobility was then remeasured, but 
likely after the native oxide had partially or fully regrown due to a time delay of several 
days. TEM analysis was carried out on several devices after transport studies were 
complete to check for the presence of stacking fault defects. Devices 1 and 3 were found 
to have zero and one fault, respectively, whereas a fourth device {D = 55 nm) with low 
mobility was found to have an atypically large fault density (see section |5] below). TEM 
analysis was not performed on device 2. 

Figure [2]^a-c) shows conductance G = Isd/Vsd, where Igd and Vsd are the source- 
drain current and bias, respectively, versus backgate voltage Vg for devices 1, 2 and 

3 at selected temperatures. The bias is set to Vgd = 1 niV. For all three devices, the 
maximum transconductance ( ) is seen to decrease as temperature is raised above 

V ^ / max 

~ 30 — 50 K. Figures [2[d-f ) show the pinchoff threshold voltages Vt corresponding to 
the data in figures |2]^a-c), where Vt is defined as the intercept between the maximum 
slope tangent hne and the G = axis. Vt typically shifts toward more positive gate 
voltages as temperature decreases, and saturates below ~ 50 K. All temperature sweeps 
reported here were from low to high temperature. We fit the pinchoff threshold data to 
an empirical function based on thermal activation Vt = Vq + Vie"'^'"''^^, where k is the 
Boltzmann constant, typically yielding an ~ 5 — 30 meV. Smaller diameter nanowires 
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Figure 3. (a) Comparison of the field-effect and effective mobilities for device 2 at 
T = 40 K. (b) The temperature dependence of effective mobility for device 2 at different 
values of gate voltage relative to Vpeak, the gate voltage at which peak mobility occurs. 
The values at Vpeak are shown by black dots, at Vpeak + 0.25 V by red dots, etc. The 
mobility at Vpeak + 0.5 V (green stars) is near the crossover point between the two 
slopes seen in the effective mobility in the left panel. 



typically have smaller Ea, such as device 3 in figure |2j^f), with Ea ^ 5 meV. Note that 
for device 1 in figure we also plot the Vt measured after the chemical treatment 
was applied to the FET channel. Vt shifted considerably to more positive gate voltage 
post-treatment, and also showed a weaker temperature dependence. This suggests that 
the surface potential and density of conduction electrons in the nanowire are controlled 
in large part by the surface chemistry [7]. 

From the measured conductance versus backgate voltage curves, both the field- 
effect mobility and the effective mobility [I3j may be extracted. The field-effect mobility 
is a lower bound on the effective mobihty, and is defined as 

1 dcr L dG . . 

where a is conductivity, n is the electron concentration, q is electron charge, C'g is the 
gate capacitance per unit length. Equation [T] only strictly holds at peak mobility, where 
= 0. The effective mobility is defined as 

LG 



where Vt is the pinchoff threshold voltage defined previously, and the expression only 
holds for Vsd « Vg — Vt. The two measures of mobilities are compared in figure [sj^a) 
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Figure 4. (a) Experimental peak effective mobilities versus temperature for devices 
1 — 3 (diameters 71, 50, and 35 nm, respectively). The empirical fitting function 
described in the text (solid lines) is given by /i = AT^{1 + Be~^'^/^'^)^'^ , where cc, 
A, B and Ea are fitting parameters given in the main text, (b) Comparison of peak 
effective mobilities versus temperature for device 1 before and after an ammonium 
sulfide etching and passivation treatment was applied to the FET channel. The fitting 
function is of the same form. For comparison, the pinchoff threshold voltages before 
and after treatment are shown in figure ^d). 



for device 2 at 40 K. The effective mobility is typically a smoother function of V^, and 
/^eff > A^fe for all of our data. Two regimes can be clearly seen in /Xefr: the slope |^^| 
is larger from Vg = —0.251^ to = +0.251^ than at more positive gate voltages. In 
figure |3j^b) we show the effective mobility versus temperature for device 2 at different 
values of gate voltage relative to the position of peak effective mobility (Vpeak)- The 
data shown are for Vg = Vpeak + ^, where the top curve (black dots) is for 6 = 0, and the 
lower curves (red, green, blue) are for 6 = 0.25, 0.5, 1.0 V, respectively. The temperature 
dependence is most pronounced at peak mobility, but follows a similar trend for points 
on the high slope region of the effective mobility curve. At large positive gate voltages 
relative to Vpeak, the mobihty shows little to no dependence on temperature. 

Devices 1 and 3 show qualitatively similar behaviour to device 2, as shown in fig- 
ure |4|^a). At a given temperature, the mobility increases with nanowire diameter, as 
was reported previously [13]. This is consistent with the mobility being dominated by 
surface charge scattering, as the overlap of the carrier distribution with the scattering 
potential becomes much stronger at smaller diameters [30]. Motivated by this hypoth- 
esis, the data in figure |4] are fit to empirical function of the form /i(T) oc T^N{T)~y, 
where N{T) is the number of surface scatterers. This function does not result from an 
analytical solution of the surface scattering problem, which is in general too difficult to 
solve without resorting to numerics [30] . Rather, this function provides a good model for 
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our data and is based on the the following reasoning. For a fixed number of scatterers, 
the average mobility increases with temperature as T^, where a; ~ 1, since the carrier 
concentration increases with temperature leading to an increase in the Fermi velocity, 
which reduces the scattering probability [30l |3T] . On the other hand, an increase in 
the number of scatterers decreases mobility. In the limit of a low density of scatterers 
and a high probability of scattering per defect, scattering events can be treated as un- 
correlated, and fi oc A^~^. However, for scattering from surface states, there is a high 
density of scatterers with a low probabihty of scattering per defect, leading to correlated 
scattering [32] (see section below). Here, the electron wavefunction remains coherent 
while interacting with multiple surface charges simultaneously, which leads roughly to 
fi oc A^~^, since the scattering matrix element is roughly proportional to A^, so the tran- 
sition rate is proportional to A^^. We model N{T) based on the thermal activation of 
surface donors: N(T) oc (1 + Be~^°-^^'^)^ where B and Ea are free parameters, similar 
to the expression used in figure [2] to model the pinchoff threshold voltages. 

The data in figure g are fit to = AT^(1 + Be''^-'^'^)-^ . For D = (71, 50, 35) nm, 
the fit parameters (excluding scaling factor A) are the following: x = (1.0,1.25,0.67), 
B = (13.4, 14.6,3.0), and Ea = (17.2, 15.1,8.0) meV. We note that the data can be fit 
equally well to a functional form /i oc A^~^, albeit with different fit parameters, but we 
chose the A^~^ form for consistency with the numerical modelling results in the next 
section. The Ea values suggest thermal ionization of the surface donor states with acti- 
vation energies in the range 8 — 20 meV, consistent with the range of Ea values obtained 
from fitting Vt in figure [2j The smaller value of B for the 35 nm diameter nanowire is 
consistent with the weaker temperature dependence of its pinchoff threshold voltage in 
figure [2]^f), indicating a smaller number of thermally activated donor states relative to 
the larger diameter nanowires. Figure |4|^b) compares the data for device 1 before and 
after an ammonium sulfide etching and passivation treatment was applied to the FET 
channel. The best fit parameters in the latter 0.62, B = 4.5, Ea = 20.2 

meV. After the chemical treatment, the turnover in mobility broadens and shifts to 
higher temperatures. This is accompanied by a much weaker change in the pinchoff 
threshold voltage with temperature, shown in figure ^d). The smaller value of fit 
parameter B after chemical treatment is consistent with the weaker temperature de- 
pendence of pinchoff threshold voltage after treatment. We note here that the detailed 
condition of the nanowire surface post-treatment is not known, and it is likely that the 
native oxide partially or fully regrew before or during the post-treatment transport mea- 
surements. The data shown here only indicate that the nanowire transport properties 
are significantly altered by chemical removal of the oxide followed by unspecified surface 
chemical processes; these processes evidently incur some change in the nature or density 
of surface states. 
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4. Numerical modelling 

We carried out numerical modelling of the nanowire transistor to test whether scattering 
from charged surface states can account for the magnitude and temperature dependence 
of the experimental mobilities. The nanowire transistor was simulated using a finite- 
element method implemented in the COMSOL® multiphysics package. The model 
consisted of a 1 fim long, 50 nm diameter nanowire atop a 175 nm thick Si02 layer with 
underlying backgate. The layer above the Si02 that embeds the nanowire is vacuum, 
with er = 1, and we take = 15.15 for the InAs nanowire. In consideration of the low 
effective mass of electrons in InAs, we used a self-consistent Poisson-Schrodinger solver 
[33] to calculate the electrostatic potential and charge distribution in the nanowire so 
that quantum confinement is properly taken into account. The model assumes that 
the conduction electron concentration at zero gate voltage is due to a surface density 
of positively charged donor states, a+ ~ 10^^ — 10^^ cm~^, an input parameter that is 
allowed to vary with temperature. 

Consider a Cartesian coordinate system with z aligned with the nanowire axis and 
radial coordinates {x,y). The potential V{x,y,z) that is a solution to the Poisson 
equation is nearly independent of axial coordinate z, so we solve the Schrodinger 
equation in a two-dimensional cross-section of the nanowire to obtain the radial 
eigenstates ipi{x,y). The electron density in the radial plane, n2d{x,y), is calculated 
from these solutions as 

POO 

n2d(x,i/) = y;n2<i,.(x,i/) = / f{E-E,)g{E-E,)dE (3) 

^0 



where g{E — Ei) = (^e-e ) one-dimensional (ID) density of states per unit 

length, f{E — Ei) is the Fermi-Dirac distribution, and Ei and il)i{x,y) are the energy 
and wavefunction of the i^^ eigenstate, respectively. The summation reduces to 

n2d{x,y) = ^v/2;;^^5^|V'^(x,i/)P x F_i/2 (^j^) (4) 

i \ B / 

where Ep is the quasi- Fermi level, Fj is the Fermi-Dirac integral of order j, and m* = 
0.023 me. The electron concentration in three dimensions is n{x,y,z) = n2d{x,y)/L, 
where L is the nanowire length. The Fermi energy Ep is determined by the net 
conduction electron concentration at zero gate voltage. Figure [sj^a) shows the values 
of o"^(T) used in the simulations, and the resulting spatial average conduction electron 
density (n(T)). We chose a function cr^(T) = ao + aie~^°-^^'^ to model the thermal 
activation of surface donor states, where ctq = 1.7 x 10^ cm~^, ai = 9.8 x 10^° cm~^ 
and Ea = 6.7 meV for the curve in figure [5|^a). These values were chosen so that the 
simulated electron density at zero gate voltage would roughly match the experimentally 
measured carrier density of device 2 at peak mobility. Note that peak mobility occurred 
at negative gate voltages in the real device, so the actual densities of surface donor states 
are likely larger than the values used in simulation. The reason for carrying out the 
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Figure 5. (a) The values of surface donor density, (T+(T), used as inputs for the 
numerical simulation of a 50 nm diameter nanowire are shown on the right vertical 
axis. The functional form, described in the text, models a simple thermal activation 
of donors. The resulting average conduction electron densities, {n{T)), are shown on 
the left axis. The cr+(r) values were chosen to produce (n(T)) at 1^ = similar 
in magnitude to the values observed experimentally for device 2 at peak mobility, 
(b) Fermi wavenumbers fci, ..,fcg of the first six radial subbands calculated from the 
Schrodinger-Poisson solutions for inputs cr+ (T). (k) is the average value over thermal 
occupation, and is proportional to the average electron velocity. 



simulations at zero gate voltage was to model the behaviour for a radially symmetric 
wavefunction, unperturbed by the presence of a nonzero gate voltage, for simplicity. 

Mobility is calculated using a multi-subband momentum relaxation time 
approximation [31]. We define three-dimensional eigenstates \m,k) = il)m{x.iy)e^^^ / y/L, 
where m is the radial subband index and k is the axial wavenumber. The transition 
probability TJ!'^ between the states |m, A;), \n,k') are calculated using Fermi's golden 
rule: 

ZjIT 

tj:;:> = j-\mj^j:,Me, - E,,) (5) 

where is the scattering matrix element (fc, m\Vc\k' ,n) resulting from the Coulomb 
interaction potential Vc of charged surface impurities. In our numerical simulations, Vc 
is obtained directly from the Poisson solver, and this takes into account both screening 
and dielectric mismatch effects [35], |36|. In the absence of these effects, Vc would be 
analytically expressed as a sum over unscreened point-charge potentials. In a cylindrical 
coordinate system (r, 6', z) where r and z are the radial and axial coordinates, 

Vc = Y. = + - cos 0. + - z,f) (6) 

where Vc.i is the potential due to a single impurity located at = {D /2,6i, Zi). With the 
numerically-derived Vc that includes screening effects, we find that the value of MJ^^, 
for a single positively charged surface impurity is on the order of 10~^ meV or less. Its 
smallness is due to the vanishing of {ipl"^ at the surface, the large dielectric constant 
for InAs, screening effects, and that the scattering potential is attractive. In this case. 
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treating scattering from single impurities independently and incoherently adding their 
rates can only lead to the observed mobilities if the surface impurity charge densities 
are unreasonably high, ~ 10^^ cm~^. At such densities, the mean separation between 
scatterers is too small for the picture of uncorrelated scattering to be valid. On the other 
hand, for a Vc that is the collective potential corresponding to a random distribution 
of many scatterers over the length of the nanowire, we are able to obtain the observed 
mobilities at impurity densities N{T) ~ cri(7') (see Figure 6). This approach justifies 
the empirical expression oc A^~^ used in the previous section to fit the experimental 
data, since the scattering matrix element M™" now roughly scales with A^, rather than 
being independent of A^ in the picture of uncorrelated single-defect scattering. 
The scattering matrix element is given by 

l>D/2 /.27r /.L/2 

Mfc"P= / / / r^,r,ir,e)Vc^Plir,e)e-'^''-'''>dzdedr (7) 

Jo Jo J-L/2 

where Vc is the total potential corresponding to a set of impurities. The integral in 
equation [7] has no straightforward analytical solution, so is generally solved numerically 
|30j . The geometry for simulating correlated scattering is indicated schematically in 
figure |6|^a), and the Poisson solution Vc obtained for a random impurity distribution is 
shown in figure |6](b). The relaxation rate in subband m due to scattering into subband 
n is calculated as 

l/r-(fc) = 5^(1- cos 0)T-^ (8) 

k' 

where is the angle of deflection between the incoming wave vector k and the outgoing 
wave vector k' . The values of k' are given by energy conservation, Em + h'^k'^ /2m* = En+ 
h'^k'^/2m* = Bp- In a ID geometry, only back scattering events contribute to electron 
relaxation rates. When the electron concentration permits the occupation of multiple 
subbands, the relaxation rate in the m}^ subband is obtained as 1/Tm{k) = J2n 
where k is the initial momentum. At low temperatures, it is valid to only consider the 
relaxation time for an electron with Fermi wavenumber kp. Making this approximation, 
we substitute the Fermi wavenumber in each subband for k. The average relaxation 
time is given by r = TiHi/n, where rij is the population of i^^ subband, leading to an 
average electron mobility jj, = er/m*. Figure [sj^b) shows the Fermi wavenumbers of the 
first few radial subbands calculated from the Schrodinger-Poisson solutions for input 
donor densities cr+(T). The first excited subband appears near 40 K, producing a dip 
in the average wavenumber (kp). The sharp drop in Fermi velocity as temperature is 
lowered below 40 K strongly increases the ionized impurity scattering rate, which causes 
a drop in mobility. 

We performed the scattering calculations in two ways: (i) calculating integrals 
MJP^, for the electron wavefunction and scattering potential over the entire length of the 
L = 1 fim nanowire, and (ii) restricting the problem to a subsection of the nanowire 
of length / < L. Method (ii) is motivated by the fact that the experimentally observed 
mobilities suggest a mean free path l^f ~ 100 — 200 nm, so that on average, we 
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expect an electron traversing the nanowire to experience several uncorrelated scattering 
events. In the latter picture, the scattering rate is calculated from the T^^, for 
the electron wavefunction restricted to a length / comparable to the mean free path, 
and the scattering rate for the entire length of nanowire is L/l times this rate. On 
the other hand, the probability for the electron to be in any one subsection is l/L, so 
these factors cancel. The only difference between the two cases is that the ID density of 
states gl, which appears in the evaluation of equation |8j is proportional to the subsection 
length. Hence, for an electron treated quantum mechanically on a length scale / (but 
classically on larger length scales), the density of states to scatter into is lower than if the 
wavefunction were spread across length L, increasing the calculated mobility. Therefore 
a factor L/l larger density of scatterers is required in calculation (ii) relative to (i) in 
order to produce the same calculated mobility. 

The results of these calculations are shown in figure [6] (d) shows the density 
of scatterers N obtained by calculations (i) and (ii) that reproduce the experimental 
mobilities. In calculation (ii), a variable subsection length / was chosen such that 
N(T) ^ o"^(T); these / values are plotted on the right axis. The calculated mobilities 
from (ii) are shown in figure 6(c) in comparison with the experimental values. A three- 
fold increase of over the range 40-150 K is able to explain the observed decrease in 
mobility with temperature for both calculation methods. Furthermore, the density of 
scatterers is nearly a perfect match to the assumed ionized surface donor density for 
method (ii). It is reasonable to expect that the increase of N with temperature results 
from the thermally activated ionization of surface donor states. Confinement also plays a 
role in this temperature dependence, since higher radial subbands contribute to a larger 
electron concentration near the surface, with a corresponding increased scattering rate. 
However, for a fixed N, this confinement effect is too small to cause a negative slope of 
the mobility-versus-temperature. We find that interband scattering plays a very limited 
role, giving at most a correction of order 10% to the scattering rates. As expected, the 
positive slope of mobility below 40 K follows the behaviour of the average Fermi velocity 
(figure [sj^b)) over the same temperature range, where only the lowest radial subband is 
occupied. Overall, the simulation results confirm that scattering from charged surface 
states at expected densities can explain the magnitude and temperature dependence of 
the experimental mobilities. 

5. Structural defects and mobility 

Finally, we studied the relationship between structure and mobility by performing post- 
measurement TEM on selected devices; this was motivated by the observation that a 
fraction of devices displayed significantly lower mobilities than were typical for a given 
nanowire diameter. A Focused Ion Beam (FIB) was used to remove devices from the 
substrate, after which they were placed on a holey carbon TEM grid for inspection. 
Indeed, it was observed that a 55 nm diameter nanowire with low mobility ~ 1,000 
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Figure 6. (a) Geometry used for calculating scattering from a random distribution 
of surface charges for a nanowire of total length L — 1 /im and diameter D —50 nm. 
The total scattering rate is obtained by calculating the scattering matrix elements 
over the entire nanowire in method (i), or by calculating the matrix elements over 
a subsection of length I and incoherently adding the rates from all L/l sections in 
method (ii). (b) Poisson potential Vc corresponding to the surface charge distribution 
in (a), projected onto a plane along the axis of the nanowire. (c) Comparison of the 
experimental mobilities (device 2) and the mobilities calculated using method (i) (the 
results using method (ii) are nearly identical), (d) The densities of surface charges 
N{T) that produce the calculated mobilities in (c) for both methods. The subsection 
lengths / used in method (ii), loosely identified with mean free path, are shown on the 
right axis. 

cm^/Vs had a high hnear density of stacking faults, at least ~(70 nm)~^ as shown in 
figure [7| In contrast, the highest mobility device we measured, device 1, had no visible 
faults along the entire channel length. Device 3 {D = 35 nm) was found to have only 
one visible fault as shown in figure [7} and better mobility than the D = 55 nm device, 
despite having a smaller diameter. The magnitude and temperature dependence of mo- 
bility appear to be greatly reduced in the D = 55 nm device due to the high density of 
stacking faults. Wurtzite InAs has a ~ 20% larger bandgap than zincblende InAs [37] , 
so that for electrons, stacking faults correspond to potential wells that may be as deep as 
~ 70 meV. Since these are planar defects, the reflection coefficient for an incoming plane 
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Figure 7. Stacking fault density and reduced mobility. Peak field-effect mobilities 
(left) and post-measurement TEM images (right) for device 3 {D =35 nm) and a low- 
mobility D =55 nm nanowire FET device. Stacking faults are indicated by the red 
arrows; at least 7 faults can be seen in the D =55 nm nanowire, compared to only one 
visible fault in the D =35 nm nanowire. The nanowires are imaged along the [2110] 
zone axis so that all planar defects will be visible. The solid lines show power law 
fits to T""-^ and T~° '^ for the 35 nm and 55 nm devices, respectively. No faults were 
observed along the entire channel for device 1 (average diameter 71 nm). 



wave can be a sizable fraction of unity. On the other hand, we cannot obtain theoretical 
mobilities as low as ~ 1, 000 cm^/Vs from a simple ID model of square well potentials 
at the linear defect density observed here. It is possible that the longer zincblende sec- 
tions may contain bound states that trap electrons [38] , leading to Coulomb scattering. 
Gap states that trap charges locally can arise at dislocations |3S], however, there are 
no mechanisms within the VLS growth method through which dislocations could form 
for the bare (111) oriented InAs nanowires studied here. A stacking fault is simply 
a rotation of the tetrahedral coordination for one monolayer, which leaves the lattice 
four-fold covalently bonded and free of distortion. Further investigation is required to 
clarify the origin of the surprisingly low mobilities seen here. 



6. Discussion 

While the data and modelling in sections 3-4 provide strong evidence for the dominant 
role of positively charged surface states as scatterers, it is also possible that negatively 
charged impurities, such as oxide charge traps |10], might play a role. Negative charges 
produce stronger scattering potentials [55], so that a relatively small number of impuri- 
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ties could limit the electron mobility. On the other hand, we observe that the pinchoff 
threshold voltage shifts to more positive values as temperature is reduced, but more 
positive gate voltages should lead to higher occupation of negative traps. Furthermore, 
if oxide charge traps limited mobility, then we would expect much higher mobilities in 
core-shell nanowires where the oxide surface is 10 — 20 nm away from the core. Some- 
what higher mobilities were observed in those nanowires [26], but only by a factor ~ 1.4 
compared to the best results with unpassivated nanowires reported elsewhere [13] and 
in the present work. Therefore we suspect that oxide charge traps do not play an im- 
portant role with regard to mobility in these nanowires, but further experiments on 
chemically and epitaxially passivated nanowires may test this hypothesis. We also have 
not considered surface roughness scattering as a dominant mechanism, though certainly 
it contributes to the scattering rate. Our reasoning is two- fold: (1) the observed tem- 
perature dependence does not appear to be consistent with surface roughness scattering, 
which should go as where is > [31], and (2) epitaxial core-shell passivation should 
smooth the roughness and lead to much higher mobilities, which does not appear to be 
the case. 



7. Conclusion 

In conclusion, we have shown evidence to support the hypothesis that ionized impurity 
scattering by charged surface states dominates the electron mobility in InAs nanowires 
across a wide range of temperatures. Transport measurements show a ubiquitous 
turnover in the temperature-dependent mobility at low temperatures. The decrease 
in mobility with temperature observed above ~ 50 K can be explained by a thermally 
activated increase in the number of scatterers. These results for pure InAs nanowires 
provide a benchmark to compare with the transport behaviour of nanowires passivated 
by chemical means or by an epitaxial shell. Additionally, by performing TEM measure- 
ments subsequent to transport measurements, we find an apparent relationship between 
moderately high stacking fault density and sharply reduced electron mobility. 
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